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ABSTRACT 

Aims. We aim to understand the nature of the absorbing neutral gas in the galaxies hosting high-redshift long-duration 
gamma-ray bursts (GRBs) and to determine their physical conditions. 

Methods. A detailed analysis of high-quality VLT/UVES spectra of the optical afterglow of GRB 050730 and other 
Swift-era, bursts is presented. 

Results. We report the detection of a significant number of previously unidentified allowed transition lines of Fe + , in- 
volving the fine structure of the ground term ( 6 D 7 / 2 , 6 D 5/ / 2 , 6 D 3/ / 2 , 6 D;l/ 2 ) and that of other excited levels ( 4 Fg/ 2 , 4 F 7/ / 2 , 
4 F5/ 2 , 4 F 3 / 2 , 4 D 7 / 2 , 4 D 5 / 2 ), from the £ a bs = 3.969, log7V(H°) = 22.10, damped Lyman-a (DLA) system located in 
the host galaxy of GRB 050730. No molecular hydrogen (H 2 ) is detected down to a molecular fraction of log / < —8.0. 
We derive accurate metal abundances for Fe + , S + , N°, Ni + , and, for the first time in this system, Si + and Ar°. The 
absorption lines are best-fit as a single narrow velocity component at z a b s = 3.96857. The time-dependent evolution 
of the observed Fe + energy-level populations is modelled by assuming the excitation mechanism is fluorescence follow- 
ing excitation by ultraviolet photons emitted by the afterglow of GRB 050730. This UV pumping model successfully 
reproduces the observations, yielding a total Fe + column density of log iV = 15.49 ± 0.03, a burst/cloud distance (de- 
fined to the near-side of the cloud) of d — 440 ± 30 pc, and a linear cloud size of I = 520± 24 ° pc. This application of 
our photo-excitation code demonstrates that burst/DLA distances can be determined without strong constraints on 
absorption-line variability provided enough energy levels are detected. From the cloud size, we infer a particle density 
of 71h ~ 5 — 15 cm -3 . 

Conclusions. We discuss these results in the context of no detections of H 2 and C I lines (with log N(C )/N(S + ) < —3) in 
a sample of seven z > 1.8 GRB host galaxies observed with VLT/UVES. We show that the lack of H 2 can be explained 
by the low metallicities, [X/H] < — 1, low depletion factors, and, at most, moderate particle densities of the systems. 
This points to a picture where GRB-DLAs typically exhibiting very high H° column densities are diffuse metal-poor 
atomic clouds with high kinetic temperatures, Tki n ii 1000 K, and large physical extents, I <; 100 pc. The properties of 
GRB-DLAs observed at high spectral resolution towards bright GRB afterglows differ markedly from the high metal 
and dust contents of GRB-DLAs observed at lower resolution. This difference likely results from the effect of a bias, 
against systems of high metallicity and/or close to the GRB, due to dust obscuration in the magnitude-limited GRB 
afterglow samples observed with high-resolution spectrographs. 

Key words, gamma rays: bursts - galaxies: abundances - galaxies: ISM - galaxies: quasars: absorption lines - cosmology: 
observations 



1. Introduction 

Thanks to prompt response, most notably on 8-10 m- 
class ground-based telescopes, high-resolution optical spec- 
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* Based on Target-Of-Opportunity observations carried out in 
service mode under progs. ID 075.A-0603, P.I. Fiore, and 075.A- 
0385, 077.D-0661, 080.D-0526, and 081.A-0856, P.I. Vreeswijk, 
with the Ultraviolet and Visual Echelle Spectrograph (UVES) 
installed at the Nasmyth-B focus of the Very Large Telescope 
(VLT), Unit 2 - Kueyen, operated by the European Southern 
Observatory (ESO) on Cerro Paranal in Chile. 



troscopy of rapidly fading long-duration 2 s) gamma- 
ray burst (GRB) afterglows provides detailed informa- 
tion on the kinematics, chemical abundances, and physi- 
cal state of the gaseous component of circumburst, host 
galaxy interstellar and i ntergalactic matter in tervening 
these lines of sight (e.g., iJakobsson et al.l [2004). At the 
GRB host- galaxy redshift, large amounts of ne utral, singly, 
and h ighly ionized gas are usually seen (e.g., iFiore et al.l 
120051 : iProchaska etafl l2007ri iFox et all l2008h . including 
the highest neutral hydrogen column densities measured to 
date at high redshift from UV absorption-line spectroscopy 
(e.g., IVreeswijk et~aT1 12004 IChen et alJl2005t iBerger et ail 



2006) . While such high H° column densities suggest an ori- 
gin in the immediate surroundings of the GRB and/or the 
inner regions of its host, the overall H° col umn density dis- 
tribu tion of the absorbers is fairly flat (Jako bsson et al.l 
2006), with values as low as ~ 10 17 up to 10 22 6 atoms 
cm" 2 (jProchaska et al.ll2008t IWatson et aLll2006D . illustrat- 
ing the diversity of environments probed by the observa- 
tions. Metallicities measured in the H° gas vary widely 
from less than l/100 th of solar to nearly solar - with 
a tendency for low m etallicities - (jFvnbo et alJ 120061 : 
iProchaska et a l. 2007b), consistent with an origin in dwarf 
and/or low-mass galaxies and the currently favoured sce- 
nario for long-duration GRBs, the collap se of a massive 
Wolf-Rayet star endowed with rotation (IWooslevI Il993t 
iMacFadven fc Wooslevlll999HWooslev fc Hegerll20(jl . 

Several nearby long-duration GRBs have now been 
shown to be associated with supernovae explosions 
dGalama et al.lll998l:ls"tanek et al.ll2003l :lHiorth et al.ll2003t 
IMalesani et all l2004tlFian et al.l l2006fr providing the most 
direct evidence that the progenitors of long-duration 
GRBs are indeed massive stars. As most star forma- 
tion occurs within molecular clouds, it is generally ex- 
pected that the latter are the birt hplaces of many GRBs . 
However, although var ious studies (jGalama & Wiiersll200il 
iReichart &: Price! l2002f ) have argued in favour of a link 
between molecular clouds and GRBs, compelling observa- 
tional evidence such as the detection of molecular hydro- 
gen (H 2 ) in GRB optical spec tra is lacking (jVreeswijk et alJ 
|2004 ITumlinson et al.l [2007h . As absorption lines related 
to ground-state or vi brationally excited H 2 energy levels 
should be detectable (jDraine &: Haol[2002T ). the absence of 
detection so far suggests that any H 2 molecules in the im- 
mediate vicinity of a GRB are dissociated by the intense 
X-ray/UV afterglow flux. However, an Hn region around 
the progenitor, or its star cluster, would also be sufficient 
to ionize H° and diss ociate H2 up to distances of 50-100 pc 
(IWhalen et al.ll2008h . 

Unlike intervening QSO absorbers and damped Lyman- 
a (DLA) systems (N(E°) > 10 20 3 atoms cm" 2 ) ob- 
served in QSO spectra, strong absorption lines involv- 
ing fine-structure and other metastable levels of ions 
such as 0°, Si + , and F e + are ubiqu i tous in GRB-DLAs 
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electrons if the particle density per unit volume is high 
enough, by indirect excitation b y IR photons, or by th e UV 
flux from the GRB afterglow. In lVreeswijk et alJ (|2007h . the 
detection of the time variability of absorption lines (see also 
iDessauges-Zavadskv et alJ l2006| ) involving the fine struc- 
ture of the ground term and other metastable energy levels 
of both Fe + and Ni + towards GRB 060418 was successfully 
modelled. This demonstrated that UV pumping followed 
by de-excitation cascades is t he mechanism at play t o form 
these lines as suggested by IProchaska et alJ (|2006l ). The 
burst/DLA distance inferred in the above case, d = 1.7 kpc, 
represented the first determination of the distance of the 
bulk of the absorbing neutral material to a GRB explosion 
site. This is consistent with lower limits on the distance, 
d 50 pc, derived from the observation of strong associ- 
ated Mgi absorption in, e.g., GRB 051111, which exclude s 
significant photo- ionization effects (jProchaska et al.ll2006h . 



Searches for molecular hydrogen in GRB-DLAs, thus 
equivalently in the interstellar medium (ISM) of GRB host 
galaxies, have led to negative results till now, though in 
two cases, towards GRB 050401 and GRB 060206, tenta- 
tive evidence of H 2 lines has been reported (IWatson et alJ 
l2006t iFvnbo et all I2006D . The H 2 column-density upper 
limits derived in a handful of cases from the non-detection 
of UV absorption lines are somewhat surprising given the 
very high H° column densities observed in some GRB- 
DLAs. Indeed, in the Galactic ISM su ch clouds are found 
to be nearly fully mol ecular (see, e.g.. ISavage et al.l 119771 : 
iJenkins &; Shavalll979l ). However, several factors could re- 
duce the amount of H 2 present in DLAs: low metallicity, low 
particle density, strong ambient UV radiation field, and/or, 
in the case of GRBs, th e UV flux from the afterglow itself. 
ITumlinson et al.l (|2007h argue that the influence of the lat- 
ter is negligible and that a combination of low metallicity 
and extreme UV radiation field from nearby star-forming 
regions - expected to be present in GRB hosts - must be 
invoked to ex plain not detecting H 2 along GRB lines of 
sight (see also lWhalen et aLll2008l) . However, given the ex- 
tremely small GRB-DLA sample size and the handful of 
systems considered so far, it is still unclear whether H 2 in 
GRB host galaxies is really deficient compared to QSO- 
DLAs and, if so, why. 

In this paper, we present an independent, thorough 
analysis of high-quality UVES spectra of GRB 050730 se- 
cured by Fiore et al. (ESO prog. ID 075.A-0603) at two con- 
secutive epochs after the burst. This line of sight exhibits 
one of the highest H° colu mn density DLAs ever observed 
in GRB afterglow spectr a ([Starling et alJl2005l IChen et ail 
120051 : iD'Elia et afll2007h with a metallic ity at the low end 
of the distribution for GRB hosts (see Chen et all [2005; 
IProchaska et al.ll2007at ID'Elia et al.ll2007F see also this pa- 
per, Sect. 1 2. 2. T|) . For the first time, we identify in the spec- 
trum of this GRB afterglow absorption lines that originate 
from numerous Fe + metastable energy levels. In addition, 
we self-consistently model the time-dependent evolution of 
the energy-level populations. This allows us to determine 
the distance of the bulk of the neutral gas to the GRB 
explosion site and the size of the absorbing cloud, and to 
address the question of whether the GRB 050730 afterglow 
flux is responsible for the observed stringent upper limit 
on the H 2 molecular fraction. We then discuss in a more 
general context the non-detections of H 2 molecules taking 
advantage of the current sample of seven VLT/UVES spec- 
tra of Swift-era, z > 1.8, GRB host absorbers. We compare 
them to QSO-DLAs where H 2 has been searched for at high 
spectral resolution to gain insights into the physical condi- 
tions prevailing in GRB-DLAs. 

This paper is organised as follows. In Sect. 12.11 we re- 
call basic facts on the observations and explain the data 
reduction process. In Sect. 12.2} we present a detailed anal- 
ysis of absorption lines from the GRB host galaxy to con- 
strain the metallicities, H 2 content, and photo-excitation of 
Fe + in the absorbing cloud. We model the time-dependent 
evolution of the Fe + metastable energy-level populations 
in Sect. 12.31 In Sect. 02 we use H° photo- ionization and H 2 
photo-dissociation calculations to assess the influence of the 
GRB afterglow radiation and the non-detections of H 2 in 
GRB-DLAs in the most up-to-date sample of VLT/UVES 
spectra. This enables us to discuss the physical nature of 
GRB-DLAs. Finally, we summarise our findings and com- 
ment on future prospects in Sect. H] 



2. GRB 050730 

2.1. Observations and data reduction 

GRB 050730 was observed wi th the Ultraviolet a nd Visual 
Echelle Spectrograph (UVES: lDekker et alJ[2000h mounted 
at the Nasmyth-B focus of the ESO VLT, UT 2 - Kueyen, 
8.2 m telescope on Cerro Paranal observatory. The data 
were gathered in service mode through a classical Target-of- 
Opportunity request by the programme of Fiore et al. (ESO 
prog. ID 075.A-0603). Two 3000 s exposures were taken 
consecutively as soon as possible after evening twilight on 
July 31, 2005, at 00:32 UT and 1:27 UT (mid-exposure 
times), corresponding to, respectively, 4.57 and 5.48 hrs af- 
ter the burst. The approximate afterglow I-band magnitude 
at the mid-exp osure time of these two epochs is I = 17.6 
and I — 18.0 (jPandev et alJ [2006) . Standard instrument 
configurations with, respectively, Dichroic #1 and Dichroic 
#2 were used (both of them using the Blue and Red spec - 
troscopic arms simultaneously; see also iD'Elia et alJ l2007) . 
The first-epoch observations cover the wavelength ranges 
of 3050 - 3850 and 4800 - 6800 A and the second epoch 
covers the 3850 - 4800 and 6800 - 10 000 A intervals. There 
are small overlaps between the two-epoch spectra and, in 
addition, two ~ 100 A gaps in the Red due to the physical 
separation between the two Red CCD detectors. CCD pix- 
els were binned 2x2, and the spectrograph entrance slit 
width was fixed to 1" to match the ambient seeing condi- 
tions leading to a resolving power of w 7 km s _1 (FWHM). 

The data were reduced using standard techniques for 
echelle data processing with the public version, v2.2.0, of 
the UVES pipeline based on the ESO MIDAS data re- 
duction package. The main characteristics of the UVES 
pipeline are to perform robust inter-order background sub- 
tractions for master flat-fields and science frames, and to 
allow for optimal extraction of the object signal performing 
sky subtraction and rejecting cosmic rays simultaneously. 
The pipeline products were checked step by step. The wave- 
length scale of the reduced spectra was then converted to 
vacuum-heliocentric values. Due to the variable nature of 
the source, individual exposures were not co-added. The 
typical signal-to-noise ratio per pixel achieved in each of 
the two reduced afterglow spectra lies in the range 8-15 
over the 4600 - 9000 A interval. 



2.2. Data analysis 

The GRB 050730 line of sight is characterized by the pres- 
ence of an extremely strong DLA absorber at z = 3.969. As 
no higher redshift system is observed and Fe II lines from 
metastable levels are also detected, we assume this is the 
GRB host-galaxy redshift. The total neutral hydrogen col- 
umn density of the GRB 050730-DLA has been measured 
in previous analys es (|Starling et alj|2005l IChen et al.l l2005: 
ID'Elia et al.ll2007t ) and the standard deviation of the three 
available measurements, based on different datasets, is very 
small (0.04 dex). We thus adopt their average value and 
combined uncertainty, i.e., logV(H°) = 22.10 ±0.10. This 
is fully consistent with our own measurement based on the 
UVES data. Numerous metal absorption lines from neutral, 
singly, and highly ionized species associated to the GRB- 
DLA are detected in the UVES spectra. While ID'Elia et all 



Table 1. Ionic column densities of neutral and singly ion- 
ized species at the GRB 050730-DLA redshift of z a b s = 
3.96857. 



Ion Energy level 


Transition 


1 „ a j i b 

log N ± (Jlog N 




lines used (A) 




1st epoch of observations: 




C° g.s. a 


1277 


<s. lO.UO 


N° g.s. a 


963, 1134.1 


14.77 ±0.03 


S+ g.s. a 


1250, 1253 


15.11 ±0.04 


Ar" g.s. a 


1066 


14.36 ±0.05 


Ni+ g.s. a 


1317 


13.53 ±0.05 



2nd epoch of observations: 
Si+ g.s. a ( 2 P 1/2 ) 1808 
Si+ 2 P 3/2 1816, 1817 

Ni+ g.s. a 1370, 1751 



15.47 ±0.03 

< 15.55 d 
13.69 ±0.02 



a ground state. 

b Listed uncertainties are the formal errors provided by 

FITLYMAN. 

c 3<r upper limit. 

d Blend with telluric absorption. 



1 http:/ /www. eso.org/projects/dfs/dfs-shared/web/vlt/vlt- 
instrument-pipelines.html 



()2007f ) have studied the multi-component structure of the 
absorption profiles from a number of transition lines us- 
ing the same dataset, we will restrict our analysis to weak 
or at most mildly saturated lines. For instance, CllA1334, 
Cn*A1335, OiA1302, Oi*A1304, Oi**A1306, SinA1304, 
and Si II* A 1309 ground-term fine-structure lines are de- 
tected but they are heavily saturated and, therefore, cannot 
be used to derive accurate column density or velocity infor- 
mation. 

The overall metal absorption profiles from neutral and 
singly ionized species at the GRB host-galaxy redshift are 
dominated by a single narrow velocity component at z a b s = 
3.96857. Interestingly, this component alone represents ^ 
90% of the total optical depth of non-saturated lines. This 
makes the Voigt-profile modelling of non-saturated lines 
from the GRB 050730-DLA cloud, as shown below, excep- 
tionally simple and robust. 

2.2.1. Metallicities 

We first focus on neutral and singly ionized species ex- 
cluding Fe + , w hich we discuss in Sect. 12.2.31 Among pre- 
vious analyses (Chen et all 120051 : iProchaska et al.l l2007al 
ID'Elia et al.l [2007) , column de nsity measurements from un- 
saturated lines were derived by IProchaska et all (|2007aD for 
S+, N°, and Ni+. In this work, we perform simultaneous 
Voigt-profile fits to the Sill, Sin*, Sn, Ni, An, and Mil 
lines detected in the UVES spectrum and, with the ex- 
ception of Sin* (see below), weak and unblended enough 
to allow for an accurate determination of column densi- 
ties. The GRB afterglow continuum was normalized locally 
around the transition lines of interest using spectral re- 
gions with a typical width of 1000 km s _1 . In order to 
fit the lines, the redshift and the (assumed-to-be) purely 
turbulent-broadening parametre values were required to be 
the same for all lines and th e two ep o chs o f observations. 
The atomic data compiled by iMortonl (|2003l ) were used for 
all studied transitio ns, except for the oscil lator strengths 
of NinAA1317,1370 ([Jenkins fe T ripd I2006D and the rest- 
frame wavelengths of S nAA1250,1253 (|Mortonlll99l . The 
spectra and best Voigt-profile fits are shown in Fig. [T] and 



0.5 - 



0.5 -I 



0.5 



0.5 




0.5 - 



0.5 - 



200 



Relative velocity (km 



Fig. 1. Normalized UVES spectral portions around weaker 
lines from neutral and singly ionized species at the 
GRB 050730-DLA redshift of z abs = 3.96857 (taken as the 
origin of the velocity scale) . Transition lines with rest-frame 
wavelengths shorter than 1350 A were observed at the first 
epoch and the other lines at the second epoch. The best 
Voigt-profile fit is overplotted. The column density from 
NinA1317, detected in the first-epoch data, has been deter- 
mined independently of the other two Ni II lines observed at 
the second epoch. The same turbulent-broadening parame- 
trc value has been used for all lines and epochs. The shaded 
regions in the panels showing the Sin lines indicate the 
expected location and shape of atmospheric features from 
a synthetic telluric absorption-line spectrum described in 
Sect. 13X51 



the results summarised in Table [JJ The best-fit turbulent- 
broadening parametre value is b — 6.5 ± 0.2 km s . The 
column density measurement uncertainties given in the ta- 
ble are the formal errors provided by FITLYMAN. They do 
not include possible additional uncertainties related to the 
continuum placement. Throughout the pa per, we adopt th e 
sol ar system abundan ces recommended in lLoddersI (|2003h . 

IChen et ahl (|2005l ) estimated the overall metallicity of 
the GRB 050730-DLA from the column density of sulphur, 
lo gjV(S+) = 15.34 ±0.09, f urther refined to be 15.22 ± 0.06 
bv lProchaska et ahl (|2007al ) . Although not considered in the 



analysis of lD'Elia et al.l (|2007l ). S II lines are also detected in 
the UVES dataset. In the latter, the 1250 A feature reaches 
zero residual intensity (see Fig. [TJ), which is not the case 
in the MIKE spectra acquired practically simultaneously. 
While the UVES spectrum is affected by a narrow spike at 
the location of the SnA1250 line, it is generally of signif- 
icantly higher quality in terms of S/N ratio and spectral 
resolution as can be s een from the compariso n of Fig. [1] 
with figs. 2 and 3 from IProchaska et alj (|2007af ) (or equiv- 
alently fig. 2 from lChen et al J 12005). A simultaneous fit of 
the S nAA1250,1253 lines in the first-epoch UVES spectrum 
leads us to derive a column density of 15.11±0.04, and a cor- 
responding metallicity of [S/H] = —2.18 ± 0.11, somewhat 
lower than the p r evious estimates. The higher estimate by 
IProchaska et al.l (|2007aD is probably caused by the use of 
the apparent optical depth method on the blended S IIA1253 
line (see Fig. QJ). 



Our result for sulphur is consistent with the Si + metal- 
licity we derive from the second-epoch UVES observations: 
[Si/H] = — 2.17±0.10. However, we note that there could be 
some contribution to the total Si + column density from the 
2 P 3 / 2 fine-structure level of the ground term. The expected 
position of the Si II* A1816 line is badly blended with a tel- 
luric absorption feature (see Fig. [TJ see also Sect. 12.231 for 
details on the telluric- line spectrum). Therefore, we shall 
consider in Table [1] the result of Voigt-profile fitting of this 
line strictly as an upper limit. Given the strength of the 
blending feature the actual Si + fine-structure level column 
density could be at least an order of magnitude smaller than 
the value reported in Table [TJ The 3cr column-density up- 
per limit derived from the non-detection of the much weaker 
A1817 line is 0.15 dex higher than this value (< 15.7). 



Singly ionized nickel is detected in the UVES spectra at 
the two epochs of observations. Ni + column densities were 
derived independently, from the 1317 A line detected in 
the first-epoch spectrum, and from the other two Ni n lines 
(AA1370,1751) observed at the second epoch (see Tabic [TJ). 
The resulting values differ by 0.16 dex (3a) in the sense 
of column densities becoming higher with time. However, 
given that the uncertainties listed in Table [TJ do not in- 
clude errors on the continuum placement, this result has 
lower significance. The column density derived from the 
Nili lin e equivalent widths measured by IProchaska et al.l 
(l2007aT ). using the sum of their three spectra taken on 
average 4.82 hrs (mid-exposure time) after the burst, is 
logiV(Ni+) = 13.56 ± 0.06. This is consistent with our 
measurements. Finally, we note that, due to the lack of 
spectral coverage of the corresponding transition lines, the 
metastable-level pop ulations of singly ionized nickel (see 
IVreeswiik et alJ[2007l ) cannot be constrained. 



Transition lines from neutral nitrogen and neutral argon 
are detected in the first-epoch UVES spectrum in a crowded 
Lyman-a forest but the lines are narrow and well-defined. 
We are therefore confident that the measured column den- 
sities (see Table [JJ are reliable. The N° column density a lso 
agrees with the one derived bv IProchaska et ahl (|2007af ) in 
the averaged MIKE spectrum. Finally, the strong C IA1277 
transition line is undetected in the UVES spectrum down 
to a 3cr column-density upper limit of 13.05 (see Fig.[JJ). 
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Fig. 2. Normalized UVES spectra around several strong 
transitions to the J = and 1 rotational levels of the vi- 
brational ground-state of H2 molecules at the GRB 050730- 
DLA redshift of z a b s = 3.96857 (taken as the origin of the 
velocity scale) in the first-epoch observations. 



2.2.2. H 2 content 

From the high quality of the UVES data, we can set a 
stringent upper limit on the presence of molecular hydro- 
gen in the GRB 050730-DLA. H2 in its vibrational ground- 
state is undetected at 3cr confidence down to column den- 
sities of logiV(H2) = 13.30 and 13.65 for, respectively, 
the J = and 1 rotational levels in the first-epoch spec- 
trum taken 4.6 hrs (mid-exposure time) after the burst (see 
Fig. [2]). The expected positions of the Lyman- and Werner- 
bands of H2 are not covered by the second-epoch observa- 
tions. We used the oscillator strengths fro m the Meudon 
groin. 3 based on calculations described in lAbgrall et al.1 
(|!994h . These upper limits translate to a molecular fraction, 
/ = 2iV(H 2 )/(2iV(H 2 ) + iV(H )) < 10- 8 , for the sum of 
the first two rotational levels. This is a fa ctor of ten deeper 
than t he constraint previously derived bv lTumlinson et alJ 
(2007}) from the aforementioned MIKE data. This is also the 
tightest constraint ever obtained for any extragalactic line 
of sight at high redshift includ i ng the database of 77 QSO- 
DLAs from iNoterdaeme et all (120081 ) . We do not find evi- 
dence in the UVES spectra of vibrationally excited H 2 lines 
as would be expected for dense cold molecula r clouds illumi- 
nated by incident GRB afterglow radiation (jDraine fc Had 
l2002t l. This is consistent with the extremely low molecular 
fraction derived above from the non-detection of transition 
lines to the H 2 vibrational ground-state. 



2.2.3. Fe + ground-state and metastable levels 

In addition to the study of metal absorption lines com- 
monly observed in QSO-DLAs (see Sect. I2.2.ip . we have 
systematically searched for all absorption features detected 
at 3a confidence in the UVES spectra and cross -checked 
their p ossible identification using the line lists from lMortonl 
(2003). To avoid false identifications due to sky absorption 
lines, we employed synthetic telluric absorption-line spec- 
tra. The latter were obtained with an IDL routine based on 
the Reference Forward Model (RFMj3, a line-by-line radia- 
tive transfer code, using the 2004 edition, vl2.0, of the high- 
resolution transmission molecular absorption (HITRAN) 
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Fig. 3. Voigt-profilc fitting to transition lines from the fine- 
structure levels of the Fe + ground-term at the GRB 050730- 
DLA redshift of z a b s = 3.96857 (taken as the origin of the 
velocity scale) in the first-epoch UVES observations. The 
lower level of the transitions, for which the column den- 
sity is determined from the fits, is indicated in each panel, 
except for the lowest Fe + energy level (ground state) corre- 
sponding to 6 D g / 2 . The shaded regions in some of the pan- 
els indicate the expected location and shape of atmospheric 
features from a synthetic telluric absorption-line spectrum 
described in Sect. 12.231 



2 http://amrel.obspm.fr/molat/ 

3 http://www.atm.ox.ac.uk/RFM/ 



database (jRothman et al.ll2005h . Featured components of 
RFM include an atmospheric profile describing mean pres- 
sure, temperature, and molecular co ncentrations for typ i- 
cally 50 atmospheric layers. Following lSmette et al.1 ((2007), 
we have calculated synthetic telluric absorption-line spec- 
tra for an amount of precipitable water vapour of 2.0 mm 
(resp. 2.4 mm) and a mean airmass of 1.3 (resp. 1.6) corre- 
sponding to the conditions prevailing during the first-epoch 
(resp. second-epoch) UVES observations of GRB 050730. 

In the line identification process, we uncovered the ori- 
gin of 18 previously unidentified narrow absorption-line fea- 
tures which are transition lines from the fine-structure lev- 
els of the ground term and other metastable energy levels 
of Fe + at the GRB host-galaxy redshift of z a b s = 3.96857. 
This is confirmed by the detection of typically two or more 
well-defined lines with consistent optical depths, and a 
cross-correlation of their measured rest-frame wavelengths 
with the NIST atomic spectra databasfl 

In order to fit the lines, we first normalized the spectra 
around each line locally in an objective manner using a cus- 
tomized routine which calculates the median value of the 
afterglow continuum after rejection of significant (> 3er) 
absorption and emission features. Note that, in spectral re- 
gions heavily affected by telluric absorption, most notably 
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Fig. 4. Voigt-profile fitting to transition lines from the fine-structure levels of the ground term and other metastable 
levels of Fe + at the GRB 050730-DLA redshift of z a b s = 3.96857 (taken as the origin of the velocity scale) in the second- 
epoch UVES observations. The lower level of the transitions, for which the column density is determined from the fits, 
is indicated in each panel, except for the lowest Fe + energy level (ground state) corresponding to 6 D 9 / 2 - The shaded 
regions in some of the panels, most notably at rest-frame wavelengths of around 1640 A, indicate the expected location 
and shape of atmospheric features from a synthetic telluric absorption- line spectrum described in Sect. [2~2~3l 



at rest-frame wavelengths of around 1640 A (see below), 
this procedure tends to underestimate the level of the con- 
tinuum. Wc assumed that the broadening parametre was 
purely turbulent. Its value as well as the redshift were as- 
sumed to be identical for all lines and the two epochs of 
observations. However, the column densities were allowed 
to be different for different energy levels and, for a given 
energy level, different from one epoch to the other. The pos- 
sibility that Fe + column densities might vary with time has 
not been considered in previous analyses of GRB 050730. 
For consistency with the time-dependent photo-excitation 
modelling of the observed Fe + level populations presented 
below, which demands the inclusion of as many energy lev- 
els as possible (see details in Sect. 12. 3p . we have consid- 



ered the more complete set of atomic parametre values 
from the Cloudy input file based on the O pacity Project 
(jVerner fc Ferlandlll996l : IVerner et al.lll999h . In view of this 
modelling, care has been exercised to derive column den- 
sity measurement errors taking into account the uncertain- 
ties in the continuum placement. These uncertainties in- 
deed dominate the formal errors provided by FITLYMAN 
in the case of weak lines. As a consequence of this, and 
contrary to Sect. 12.2. 1| errors have been estimated from 
the line fitting using the best-fit continuum normalisation 
(see above) comparing the results with Voigt-profile fits 
with the continuum placed at reasonable upper and lower 
boundaries. We found these to correspond to adjusting 
the continuum level by plus or minus 0.5 times the noise 



Table 2. Column densities of the fine-structure levels of 
the ground term and those of metastable levels of Fe + at 
the GRB 050730-DLA redshift of z abs = 3.96857. 



Ion Energy level Transition log TV ± <7i og jv b 



lines used (A) 


1st epoch of observations: 






Fe+ 


p. a ( 6 D 9/2 ) 


1083, 1112 


I5.3il° ° 4 


Fe+ 


D7/2 


1266, 1267 


14.391°;!° 


Fe+ 


6 D5/2 


1271 




14.361°;?? 


Fe+ 


D3/2 


1275.7, 1275.8 


14 04 +0 12 


Fe+ 


6 D 1/2 


1277.64, 1277.68 


< 13.60 c 


2nd epoch of observations: 






re 


g.s. ( D 9/2 ) 


1611 




1 5 q 4 +0.06 


Fe+ 


6 D7/2 


1618, 


1621 


14.4ll°;°2 


Fe+ 


6° 5/2 


1629, 


1631 


13-961H? 


Fe+ 


D3/2 


1634, 


1636 


1 3 7Q+ 09 


Fe+ 


6 D 1/2 


1639 




< 13.65 d 


Fe+ 


4 F 9/2 


1559, 


1566, 1569 






1610, 
1702 


1612, 1658 


. 1 j _i_n riQ 

i4.3ilg:jg 


Fe+ 


4 F 7/2 


1563, 


1625, 1659 






1712 




13.641°;!° 


Fe+ 


4 F 5/2 


1570, 


1584, 1633 


13.501°;!! 


Fe+ 


4 F 3/2 


1574, 


1640 


< 13.39+° !^ d 


Fe+ 


D 7 / 2 


1635 




■I 4fi +0.08 
1_. -O_ . 8 


Fe+ 


4 D 5/2 


1641 




< 13.291°;!^ d 



a ground state. 

b Listed uncertainties were derived from the line fitting using 
the best-fit continuum normalisation plus or minus 0.5 times 
the noise RMS in the adjacent continuum. 
c 3a upper limit. 

d Possible blend or mis-identification (see text). 



RMS in the continuum adjacent to the line. The results 
of Voigt-profile fitting together with the list of lines used 
for the fits are shown in Figs. [3] and 31 and summarised in 
Tabled The best- fit turbulent-broadening parametre value 
is b = 9.9 ± 1.1 km s _1 . Note that the Fen lines were fitted 
separately from the other metal lines resulting in a some- 
what different b value (see Sect. I2.2TTT) . However, since the 
overall fit is based on weak lines the column density results 
are not very sensitive to the exact b value. 

Our measurements of column densities for the fine- 
structure levels of the Fe + ground-ter m ( 6 D) are i n most 
case s consistent with those d erived bv lD'Elia et al.l ([2007) 
and iProchaska et alj (|2007ah . This comparison takes into 
account the fact that the oscillator strengths we have 
adopted differ slightly (by 10-30%) from those compiled by 
iMortonl (|2003[ ). which were adopted in previous analyses. 
The only exception to this consistency is the level, 
whose A1639 transition line is not detecte d in the second- 
epoch U VES spectrum (see Fig. H]) while IProchaska et al.l 
(|2007al ) report a column density of logiV = 13.65 ± 0.07. 
This is 0.65 dex higher than our 3cr detection limit, log A < 
13.0. However, the latter may be biased low due to our 
continuum normalisation being affected by the presence of 
telluric absorption features, which are conspicuous in this 
part of the spectrum as seen from the telluric absorption- 
line template (Fig. H|) . The template also matches the four 
absorption-line features at -120, -80, 0, and +80 km s _1 
observed around 8145.5 A in the MIKE spectrum (see bot- 
tom right panel of fig. 3 in IProchaska et al1l2007ah . both in 



terms of wavelength positions and relative optical depths, 
pointing to significantly higher amount of precipitable wa- 
ter vapour present, during the observations of GRB 050730, 
in the atmosphere above Las Campanas observatory than 
on Paranal. From this, we conclude that the actual Fe + 
6 D 1 / 2 column density must be smaller than 13.65, which 
we conservatively adopt as an upper limit in Table [2] 

Taken altogether, the column densities measured for the 
fine-structure levels of the Fe + ground-term ( 6 D7/ 2 , 6 D5/2, 
and 6 D 3 / 2 ) are consistent with the values becoming smaller 
with time. However, the statistical significance of this result 
is at most marginal when levels are considered individually. 
Uncertainties are large for the first epoch of observations 
because, at A ros t < 1300 A, only a few lines are available 
for a given energy level, i.e., either a single line or two mu- 
tually blended lines (see Fig. [3]). It should thus be kept in 
mind that for the first epoch some of the fitted lines might 
be blended or simply mis-identified. In addition, these lines 
are weak and noisy. The 6 D 1 / 2 level is undetected in the 
first-epoch UVES spectrum with a 3cr column-density up- 
per limit of 13.60 (see Fig. [3] and Table [5]). In contrast, 
the second-epoch UVES spectrum displays a wealth of Fe 11 
lines from metastable energy levels. Column densities for 
the second epoch of observations are in most cases secured 
from at least two well-defined lines, up to seven in the case 
of the 4 Fg/ 2 level. Exceptions, for which we conservatively 
quote upper limits or equal values in Table are the F 3 / 2 
level (for which the feature observed at A ros t ~ 1574 A is 
inconsistent with the non-detection of the A1640 line prob- 
ably because of improper continuum placement) and the 
4 D 5 / 2 level (whose detection through the A1641 line is un- 
certain due to the presence of adjacent telluric absorption; 
see Fig. HJ). 

2.3. Fe + photo-excitation modelling 

Our group successfully modelled the time-variable exci- 
tation of Fe + and Ni + at the host-galaxy redshift of 
GRB 060418 as due to pumping by afterglow ultra-violet 
(UV) photons (jVreeswiik et al.ll2007h . This led to the first 
determination of the distance of the bulk of the host- 
galaxy absorbing neutral material to a GRB (d = 1.7 ± 0.2 
kpc). The detection of absorption- line variability with a 
steadily rising population of Fe + in the 4 F g / 2 energy level, 
and the simultaneous decay of the other levels, ruled out 
collisions with electrons as the dominant source of ex- 
citation. There is now mounting evidence that excita- 
tion by UV photons and subsequent de-excitation cas- 
cades are the mechanisms at play to populate the Fe + 
metastable energy levels in GRB afterglow spectra, as vari- 
ations of Fen absorption lines have been observed at the 
host-galaxy redshifts of two additional bursts, GRB 020813 
and GRB 08 0319B, for which suitable multi-epoch data 
are a vailable (jDessauges-Zavadskv et al.ll2006l : iLTElia et al.l 
120091 ). In the case of GRB 050730, absorption-line variabil- 
ity of transitions involving the Fe + 6 D g / 2 , & ^7/2 > 6 Ds/2 > and 
6 D 3 / 2 levels cannot be firmly established (see Sect. I2.2.3|) 
as the difference between the two epochs of observations 
is small relative to the time elapsed since the burst and, 
therefore, collisions cannot be excluded as a possible source 
of excitation. Single-epoch Boltzmann fits to both the ob- 
served ground-state and the metastable level populations 
of Fe + , including 4 F g / 2 and higher energy levels, give poor 
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Fig. 5. Top panel: Time-dependent evolution of the column densities of the ground and several metastable levels of 
Fe + in our best-fit UV pumping model (indicated with solid lines) using Cloudy atomic parametres and allowing for the 
absorbing cloud to be extended (see discussion in text). Measurements (resp. upper limits) are indicated with open circles 
(resp. triangles); the upper limits were not used to constrain the fit. Bottom panels: Zoom-in around the two epochs 
of observations for three different model fits. The bottom right panel shows the same model as the top panel but also 
identifies the Fe + levels of the observed and modelled populations. The other two bottom panels show fits to the data 
assuming a zero cloud size, where the old atomic parametre fit is featured in the left panel and the Cloudy parametre fit 
is located in the middle. 



results, with X l = 13.2/(4 - 2) = 6.6, T cx - 400 K, 
and logiV(Fe + ) = 15.3 for the first epoch, and xt = 
106.8/(8-2) = 17.8, T cx - 2600 K, and log7V(Fe+) = 14.5 
for the second epoch (while logiV(Fe + ) = 15.34±0.06 is ob- 
served). The very poor Boltzmann-distribution fits, as well 
as the disparate excitation temperatures and ground-state 
column densities inferred from the two quasi-simultaneous 
epoch data, lead us to reject the hypothesis that collisions 
are the cause of the excitation of Fe + . As a consequence, 
we apply below our UV pumping model to the case of 
GRB 050730. 



2.3.1. Model 

In our model, GRB afterglow photons progressively ex- 
cite Fe + ions in a gas cloud located at a distance d from 
the burst. Although we take into account both direct ex- 



citation by afterglow IR photons and indirect excitation 
by UV photons, the latter is, at a given distance, by far 
the dominant excitati on mechanism ( Prochaska et al.f 2006; 
IVreeswijk et al.l [20071 ). In addition, as a result of the large 
separation of the Fe + ground-term fine-structure energy 
levels, the cosmic microwave background radiation is al- 
ways a negligible sou rce of excitation of Fe + even at z = 4 
(|Silva fc Viegadl2002D . The ion level populations at a given 
time are a function of the strength and spectral slope of 
the afterglow UV radiation and its decay with time, the 
pre-burst Fe + column density (all ions are assumed to be 
initially in the ground state) , the distance of the absorbing 
cloud to the GRB explosion site, and the velocity broaden- 
ing of t he ions. In this r e spect , our analysis closely follows 
that of IVreeswijk et al.l (|2007h who modelled the case of 
GRB 060418. Although the narrow time window covered 
by the UVES observations of GRB 050730 does not allow 



for a definite detection of absorption-line variability (see 
above), unlike the UVES observations of GRB 060418, the 
large number of Fe + metastable energy levels detected in 
absorption towards GRB 050730 constrains the model quite 
well, as shown below, so that the burst/DLA distance can 
be reliably determined. 

The intensity of the UV flux illuminating the cloud is 
obtained by converting the observed afterglow light curve 
to the GRB rest frame . We used the light-curve description 
of IPandev et a l. (2006), namely a broken power law with a 
break time around 0.1 day with pre- and post-break power- 
law indices of a x = -0.60 ± 0.07 and a 2 = -1.71 ± 0.06, 
respectively. These authors also measured a spectral slope 
of = —0 - 56 ± 0.06, which we adopted in our model. 
iKann et al.l (|2007h measured similar values for these in- 
dices. For the flux zero-point of the light cur ve, we adopted 
I = 17.22 at i bs=3.47 hrs after the burst (|Pandev et al.l 
l2006| k this epoch was also assumed to be the time of the 
jet break. We note that the /-band filter does not contain 
the Lya forest of absorption lines whi ch in turn affects th e 
B, V, and i?-band photometry (see IPandev et al.l l2006h . 
while the rest-frame central wavelength of the /-band fil- 
ter is around 1600 A, which is in the wavelength regime 
(912-2600 A) of the UV photons responsible for the excita- 
tion of Fe + . This leads to the following description of the 
rest-frame flux at a distance d from GRB 050730: 
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in erg s _1 cm -2 Hz -1 . In the conversion to the rest frame, 
we a dopted the Gal a ctic e xtinction of E(B — V) = 0.046 
from Schlcgc f et alJ {1998) and assumed that any other 
extinction along th e line of sight is negligible (see also 
IStarling et al.l [2005T) . The latter assumption is consistent 
with the low metallicity and low depletion factor of [S/Fe] 
measured at the host-galaxy absorber redshift, which indi- 
cates a low dust content if any (see Sect. 12. 2 Jj) . For the cal- 
culation of the luminosity distance, di — 3.55 x 10 10 pc, we 
adopted H = 70 km s" 1 Mpc" 1 , n M = 0.3, and ft A = 0.7. 
For the calculation of the atom level populations, we re- 
fer the reader to equations 3- 5 and the accompanying ex- 
planations of IVreeswijk et all (|2007| ). We checked for the 
contribution of the source function by performing two runs 
including and excluding it; the change in the resulting chi- 
square was found to be negligible. As including the source 
function requires much more CPU time, we set the source 
function to zero. For the velocity broadening of the ions, we 
adopted the value that we measured in the UVES spectra, 
b = 10 km s _1 (see Table [5} as it is well-constrained due 
to the simplicity of the absorption-line profile of this GRB- 
DLA. For many UV transitions, th e cloud that we model i s 
optically thick and, therefore, as in IVreeswiik et~aTI (|2007l) . 
we sliced up the cloud in a sufficient number of plane- 
parallel layers, so that each layer can be considered optically 
thin for a given transition. In contrast to GRB 060418, we 
fixed the starting time to to 30 s in the rest frame for GRB 
050730 as the results are very insensitive to any value of to 
between five and 200 s. In other words, the level popula- 
tions at roughly an hour (in the rest frame) after the burst 



do not depend on the brightness of the very early afterglow, 
up to a rest-frame time of about five minutes. 

2.3.2. Atomic data 

We considered two different sets of Fe + atomic parametres 
in our calculations : the first one is the same set as used 
in IVreeswijk et aD (|2007l ). and the second one is the atom 
model employ ed in the Clo udy photo-ion i zation code (see 
iFerlandl [20031 ) described in IVerner et"aH (|l999h . The for- 
mer set of atomic data, which we refer to below as the 
"old atomic parametres" , includes the 20 lower energy lev- 
els of Fe+ (up to E = 18,886.78 cm" 1 ). The probabili- 
ties for spontaneous decay, or A values, of the forbidden 
transitions b e tween all these lower levels were taken from 
iQuinet et al.l (|1996[) . For the allowed transitions between 
lower and higher (starting from E = 38,458.99 c m" 1 ) ex- 
cited levels, we adopted the A values compiled bv iMortonl 
(2003]) w henever availa ble and, if not, we used those pro- 
vided bv lKuruczl (|2003r FI. The Cloudy atomic data set in- 
clud es all 63 lower ener gy levels of Fe + . We here directly 
cite IVerner et al.1 (|1999T) for the description of the corre- 
sponding sources: "Transiti on prob a bilitie s are taken from 
theoretical calculations bv iNaharl (|l995l allowed transi- 
tions) and IQuinet et alJ (|l996l forbidden trans itions) and 
supple ment ed by data from compilat i ons by iFuhr et ail 
(|l988f ) and iGiridhar k. Arellano Ferrol (|l995l ). Transition 
probabilities for all intercombinatio n lines not covered by 
these compilations are taken from the lKurucz fc B ell (1995) 
database. Uncertainties are generally smaller than 20% for 
strong permitted lines but can be larger than 50% for weak 
permitted and intercombination lines. For the forbidden 
transitions of interest, uncertainties are expected to be less 
than 50%." Using these two different atomic data sets gives 
us insight into the uncertainties in the output fit parame- 
tres, most importantly the distance of the photo-excited 
cloud to the GRB. The final resulting uncertainties of the 
model parametres only reflect the uncertainties in the col- 
umn density measurements and atomic data used, and do 
not include the fact that our model of a plane-parallel se- 
ries of absorption slabs is a simplification of the true cloud 
structure in the host-galaxy ISM. We also note that the dif- 
ferent spontaneous decay coefficients A u \ (or equivalcntly /) 
of each data set result in slightly different column densities 
measured in the spectra. 

2.3.3. Results 

We first ran our photo-excitation code using the old atomic 
parametres, which results in a poor fit of the data with a 
reduced chi-square of \l = 236.8/ (13 — 3) = 23.7 (see bot- 
tom left panel of Fig. [5]). The corresponding best-fit values 
for the distance of the cloud to the GRB explosion site 
and the total Fe + column density are d = 550 ± 40 pc 
and logiV(Fe + ) = 15.25 ± 0.02, respectively. When em- 
ploying the Cloudy atomic parametres, the fit greatly im- 
proves (see bottom middle panel of Fig. [5]), with xl = 
24.8/(13 - 3) = 2.5, and best-fit values of d = 470 ± 30 pc 
and log V(Fe + ) = 15.52±0.03. The la errors reported here 
correspond to the parametre value for which the overall chi- 
square increases by one with respect to the minimum. This 
parametre value is searched for by performing several new 

5 http://kurucz.harvard.edu 



fits where the relevant parametre is fixed, while the other 
parametres are left free to vary, until a fit is found where 
Ax 2 equals one. 

The difference between the best fits, obtained using 
either the old atomic parametres or the Cloudy atomic 
parametres, mainly comes from the very different values 
of the spontaneous decay coefficients of the Fe + 4 F 9 / 2 and 
4 F7/ 2 energy levels. The lower A- values in the old atomic 
parametre set, by about a factor of three, cause these lev- 
els to be severely overestimated in the model, as evident 
from the bottom left panel of Fig. [5l This in turn results 
in the populations of the ground state and associated fine- 
structure levels of Fe + being underestimated. Instead, the 
Cloudy atomic parametre set results in a very reasonable 
fit as shown in the bottom middle panel of Fig. 03 However, 
we note that despite the large difference in the atomic data 
parametre sets used (both in the number of levels and A- 
values), the derived distances: d — 550 ± 40 pc (old set) 
and d = 470 ±30 pc (Cloudy set) are consistent at the 1.6c 
level. 

In contrast to the case of GRB 050730, the observa- 
tions of GRB 060418 w ere successfully fit u sing the old 
atomic parametres fsee IVreeswiik et al.ll2007fj . Therefore, 
we verified whether the time-variability of the excited Fe + 
and Ni + energy-level populations observed along the line 
of sight to GRB 060418 can also be fit using the Cloudy 
atomic parametres. We find that the latter actually pro- 
vides a much better fit, with a reduced chi-square (xl = 
19.3/(36 — 5) = 0.62) almost half of that of the fit using the 
old atomic parametres. The Cloudy best-fit distance of the 
absorbing cloud to the GRB (d = 2.0 ±0.3 kpc) is fully con- 
sistent within the errors with the previously derived value 
(d = 1.7 ±0.2 kpc), while the best-fit total Fe + column den- 
sity is significantly higher: log7V(Fc + ) = 15.36 ±0.11 (to 
be compared to log7V(Fe+) = 14.75±^ from the old fit). 
This updated Fe + column density is consistent with that 
observed along the line of sight to GRB 060418, thereby 
eliminating the need for another absorbing cloud, located 
at much larger distance from the b urst, that is not signif i- 
cantly excited (see the discussion in IVreeswiik et al] [2007'). 
This major improvement of the photo-excitation modelling 
of the GRB 060418 data lends strong additional support for 
adopting the Cloudy atomic parametres in our code rather 
than the old atomic parametres. 

Up to this point, we have assumed that all the absorb- 
ing neutral material is located at the same distance along 
the line of sight to the GRB, i.e., that the cloud is infinitely 
small. In an improved version of our code, we allow the ab- 
sorbing cloud to have a physical extent over which the cloud 
is assumed to have constant density, and consider its linear 
size / as an additional fit parametre. Applying this more re- 
alistic model to the data, with each cloud layer (see above) 
now located at a different distance from the burst, leads to 
the following results: a burst/cloud distance (defined to the 
near-side of the cloud) of d — 440±30 pc, a linear cloud size 
(i.e., along the line of sight) of I = 5201 24 >q pc, and a total 
Fe+ column density of logiV(Fe+) = 15.49 ± 0.03. It turns 
out that the goodness-of-fit improves significantly when in- 
troducing the cloud size, with xl = 12.2/(13 - 4) = 1.4 
(compared to xl = 2.5 before), and the best-fit parametre 
values are well constrained. We will therefore adopt these 
model fit results that include a cloud size in the rest of the 
analysis. Although in reality the structure of ISM clouds 



could be much more complex than our single homogeneous 
cloud model, the simple single-component absorption-line 
profile in the case of GRB 050730 (see Figs. [TJ |1 and g]) 
indicates that our model may actually be a reasonable ap- 
proximation for this burst. The fit including a cloud size is 
shown in both the top panel and the bottom right panel of 
Fig. [5] We also note that the change in light-curve decay 
index around the break time of 3.47 hrs is evident in the 
evolution of the modelled Fe + level populations. 

3. Molecular hydrogen in GRB-DLAs 

3.1. GRB-DLA gas distance and H2 photo-dissociation 

In the previous sectio n, we determined t he second precise 
distance (but see also ID'Elia et al.ll2009f ) from a GRB to 
the bulk of the absorbing neutral material experiencing 
photo-excitation of its Fe + component by the time- variable 
GRB afterglow radiation. We found that the GRB 050730- 
DLA is located 0.5 kpc away from the GRB explosion 
site. This distance is large enough so that the gas is not 
strongly ionized by the incident radiation and, therefore, 
remains essentially neutral, in agreement with the ob- 
servations (see below). This is yet closer than what we 
previously determined towards GRB 060418 at z = 1.49, 
where the burst / PL A distance was fou nd to be 2 kpc 
(see Sect. HO] and IVreeswiik et all l2007fi . More recently, 
ID'Elia et al.l <|2009h determined distances of 2-6 kpc for 
the different absorption systems in the host galaxy of 
GRB 080 319B, using the same ty pe of modelling as we pre- 
sented in IVreeswiik et aD (2007) but with different atomic 
data parametres. Although it is an open question why 
the bulk of the absorbing neutral material is located at 
such distances from the GRB explosion sites, these results 
could be a consequence of the host galaxies of high-redshift 
long-duration GRBs being compact, possibly H°-rich, dwarf 
galaxies. This would be consistent with the low metal- 
licity of the GRB 050730-DLA which we measured to be 
[S/H] = -2.2 (see Sect. |2~2~1"T) . This would also be in line 
with the findings of imaging and spectroscopic studies that 
the typical long-duration GRB host at intermediate red- 
shift is a blue, sub-luminous, low-mass, st ar-forming galaxy 
dLe Floc'h et all 120031: IChristensen et al.l 12004 IChen et all 
120091 : ISavaglio et al.ll2009D . In the following, we assess m 
more detail the influence of the GRB afterglow radiation 
on the neutral (H°) and molecular (H2) contents of a DLA. 

Using the observational constraints on the light curve 
described in Sect. 12. '6\ we can calculate the total number of 
H°-ionizing photons in units of s _1 cm -2 at a given distance 
from the burst as: 

$™F v <T v dv/{hv) 

"-phot = 75o -7—, (2) 

\ va a v dvjv 

where F v is the rest-frame GRB afterglow flux given in 
Eq. [TJ vq the frequency of the hydrogen ionisation thresh- 
old (13.6 eV), and a v the H° photo-ionisation cross-section 
For th e latter, we adopt the fits presented in IVerner et al.l 
(|l996l ). riphot can be integrated over time from the ob- 
served light curve, providing the total column density of 
H°-ionizing photons at a distance d released by the GRB. 
The left panel of Fig. [6] (dashed line) shows the result of 
this calculation for GRB 050730 photons emitted between 
30 and 3310 s after the burst (rest frame), i.e., up to the 
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Fig. 6. Left panel: Column densities of H°-ionizing and ^-dissociating photons released by GRB 050730 between 30 and 
3310 s after the burst (rest frame; i.e., up to the mid-exposure time of the first-epoch UVES spectrum) as a function of 
distance to the burst. Right panel: H° photo- ionization and H2 photo-dissociation rates integrated over the rest-frame 
time interval 30-3310 s as a function of distance to GRB 050730, neglecting self-shielding and assuming no dust (see text). 
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Fig. 7. Neutral hydrogen (H° ; left panel) and molecular hydrogen (H2 ; right panel) column density remaining at the time 
of the first-epoch UVES spectrum as a function of distance to GRB 050730 for dust-free clouds with different initial (pre- 
burst) column densities. The initial column density (10 14 -10 23 ) of each calculation is indicated with a dotted horizontal 
line. We performed two sets of calculations: one where the cloud size is negligible (solid lines), and one where the initial 
cloud size is assumed to be equal to the burst/cloud distance (dashed lines). In each panel, the diamond symbol indicates 
our observational constraint, i.e., log7V(H°) = 22.10 on the left and logiV(H2) < 13.8 on the right, at a distance of 
d ~ 0.5 kpc from GRB 050730 as determined from the UV pumping model. 



time the first spectrum was taken. In the right panel of 
Fig. El (dashed line), we show the H° photo-ionization rate, 
or i?i on . It is calculated in the same way as above except 
that only the numerator of Eq.[2J and not its full expression, 
is considered: i?i on = / F v a v dv/(hv). 

In addition, we calculate the total number of H2- 
dissociating photons, i.e., the number of photons in 
the energy range 12.24-13.51 eV (corresponding to the 
Lyman- Werner absorption bands), that are released by 
GRB 050730 between 30 and 3310 s after the burst (rest 
frame). The result is shown as a function of distance to the 
burst in the left panel of Fig. [6] (solid line). Only 10-15% 
of these p hotons will actually lead to dissociation of H2 
molecules (jDraine fc Bertoldil [l996t iHollenbach fc Tielensl 



I l999t). The H9 p hoto-dissociation rate is given by (see, e.g., 
kbel et al.l[l997h : 

-Rdiss « 1.1 X 10 8 JlW S'shicld S 1 (3) 

where Jlw (in units of erg s _1 cm -2 Hz -1 ) is the rest- 
frame GRB afterglow flux at mean energy hv — 12.87 eV, 
and S'shicid a correction factor for H2 sel f-shielding and de- 
struc tion of UV photons by dust grains (jDraine fc Bertoldl 
199G: see also iHirashita fc Ferrarall2005l ) . In the right panel 
of Fig. [6] (solid line), we plot i?diss, or equivalently &27 
(jAbel et ~alJll997ft . integrated over the rest-frame time in- 
terval 30-3310 s as a function of distance to the burst. We 
neglect self-shielding and assume no dust as the most con- 
servative approach to reveal the maximum influence of the 
propagating GRB afterglow radiation. It can be seen from 



Table 3. Neutral gas, metals, dust, and molecules in the VLT/UVES sample of Swift-era, z a b s > 1.8 GRB hosts (as of 
June 2008). 



GRB 


Zabs 


logiV(H u ) 


[X/H] 


X 


[X/Fcf 


logiV(Fe)dust 


log A^Ha) 1 


lo, 


g/ h 


logiV(C u ) 


050730 


3.969 


22.10 ±0.10 


-2.18 ±0.11 


S 


-0.06±0.06 d 


< 14.77 d 


< 13.8 s 


< 


-8.0 


< 13.05 g 


050820 


2.615 


21.05 ±0.10 


-0.39 ±0.10 


Zu 


±0.80 ± 0.03 


16.05 


< 14.1 g 


< 


-6.7 


< 12.48 g 


050922C 


2.200 


21.55 ±0.10 


-1.82 ±0.11 


s 


±0.76 ± 0.05 c 


< 15.12° 


< 14.6 g 


< 


-6.7 


< 12.45 g 


060607 


3.075 


16.95 ± 0.03 










< 13. 5 g 


< 


-3.1 


< 12.03 s 


071031 


2.692 


22.15 ±0.05 


-1.73 ±0.05 


Zu 


±0.04 ±0.02 


14.83 


< 14.1 g 


< 


-7.8 


< 12.80 g 


080310 


2.427 


18.70 ±0.10 


< -1.91 ± 0.13 b 


o b 


a 


a 


< 14.3 g 


< 


-4.2 


< 12.55 s 


080413A 


2.435 


21.85 ±0.15 


-1.60 ±0.16 


Zn 


±0.13 ±0.07 


15.13 


< 15.8 s 


< 


-5.7 


< 13.28 s 



References to the ob servation s: 050730: ID'Elia et all l|2007f ) : 050820: iLedoux et all (|2005l ): 050 922C: iD'Elia et all ( 2005): 060607: 
iLedoux et all (2006); 071031: Eedoux et all (|2007T >; 080310: IVreeswiik et all i|2008ah ; 080413A: IVreeswiik et all l|2008bl ). 
a The metallicity and/or depletion factor cannot be determined due to ionization effects. 

b The OlA1302 line used for the fit is partially blended with strong AlllA1670 absorption from an intervening system at 
Zabs = 1.671. 

c The total Fe + column density includes the contributions of energy levels above the ground state up to 4 F 9/ / 2 whenever 
transition lines from these levels were detected. 

d When [X/Fe] is negative, N(Fe)a us t cannot be calculated directly. In this case, we estimate a 3cr upper limit on -/V(Fe)d us t by 
considering the upper bound provided by the 3cr error on [X/Fe]. 

Significant overabundance of sulphur compared to iron might be present in this system (see Sect. 13.2^ 1 . The column density of 
Fe in dust derived from [S/Fe] (as all measurements based on S or O) should thus be considered strictly as an upper limit. 
f Sum of the contributions of the J = and 1 rotational levels. 
s 3cr upper limit. 

h / = 2Ar(H 2 )/(2Ar(H 2 ) + A r (H )). 



this figure that i?diss = 1 at d « 500 pc and, therefore, that 
H2 photo-dissociation can only be effective at d < 500 pc. 

From the definition of the rate constant, i.e., dN — 
RNdt, we now calculate the H° and H2 column densities re- 
maining at the time the first GRB 050730 UVES spectrum 
was taken, i.e., 3310 s after the onset of GRB 050730 (rest 
frame), for clouds located at given distances to the burst 
and having different initial (pre-burst) column densities 
(see Fig. [7]) . Each calculation starts out by placing a cloud 
with an assumed initial column density at large distance, 
where H° ionization and H2 dissociation are negligible. As 
the distance decreases, the cloud is increasingly affected by 
the afterglow photons, until the distance is reached where 
the entire column density is ionized or dissociated. This is 
shown in Fig. [7] for both H° ionization (left panel) and H2 
dissociation (right panel). The two different line types de- 
picted in Fig. [7] correspond to two different assumptions 
for the cloud size: one assumes that the cloud size is negli- 
gible (solid lines), while the other assumes that the initial 
cloud size is equal to the burst/cloud distance. The calcula- 
tions take into account the fact that the column density of 
photo-ionized (or photo-dissociated) particles never exceeds 
the column density of photons available at a given distance. 
We here again conservatively neglect self-shielding and, in 
addition, assume no dust is present, which is relevant to 
the case of GRB 050730. 

One can see in Fig. [7J that, at typical distances of 
0.5 kpc and larger, incident GRB afterglow radiation can- 
not significantly ionize H° (left panel) nor dissociate H2 
(right panel) within a DLA. For a DLA cloud, i.e., with 
log AT(H°) > 20.3, significant ionization effects can only oc- 
cur at distances smaller than 50 pc. Moreover, at distances 
larger than 0.5 kpc the decrease in H2 column density is 
smaller than 0.4 dex. Including the process of destruction 
of UV photons by dust grains would decrease the influence 
of the GRB even more. However, this will be important 
only for high A7(H°) values coupled to high metallicities 



and large dust contents, while for GRB 05073 a low dust 
content was inferred (see IStarling et al1l2005h . Therefore, 
the current lack of H2 detections in GRB-DLAs located 
at least 0.5 kpc away from the GRB explosion sites must 
be caused by one or several factor s other than the propa - 
gating afterglow radiation (see also lTumlinson et al.ll2007l ). 
This is what we investigate below by considering the most 
up-to-date sample of high-resolution VLT/UVES spectra of 
high-redshift GRB afterglows. 

3.2. Statistics of H 2 in GRB-DLAs 

In Table [3J we list the seven Swift-era. GRB afterglows 
with redshifts higher than 1.8 for which UVES observa- 
tions were secured as of June 2008. For each line of sight, 
we measured in this work total H° column densities from 
Voigt-profilc fitting to the Lycv and/or Ly/3 absorption lines 
at the GRB host-galaxy redshifts (denoted z a b s in the ta- 
ble). In the same way as described in Sect. 12.2. 11 we have 
also measured or re-measured in an homogeneous manner, 
from Voigt-profile fitting to the associated metal absorption 
lines, metallicities, [X/H] = log JV(X)/iV(H) - log(X/H) , 
and depletion factors, [X/Fe], where X is a non-refractory 
reference element. X is taken to be Zn when Znii lines are 
detected, or S otherwise. For GRB hosts that do not meet 
the requirement for DLA absorbers (logiV(H°) < 20.3), the 
abundance of neutral oxygen may be used instead. This 
atom is linked to neutral hydrogen by a strong charge- 
exchan ge reaction due to the similarity of ionization po- 
tential ((Field fe Steigma"n1ll97lD . and the [0°/H°] r atio can 
lead to an estimate of the total oxygen abundance (|Viegasl 
Il995f ). The ionization correction factor is smaller than 10% 
for GRB 080310 having a logarithmic H° column density 
of 18.8. However, for GRB 060607 with log AT(H°) = 16.95 
the ionization correction highly depends on the hydrogen 
particle density and, therefore, no robust metallicity esti- 
mate can be derived. To measure depletion factors and get 
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Fig. 8. Le/i panel: Logarithm of the neutral hydrogen column density versus metallicity, [X/H] = log N(X)/N(R) — 
log(X/H) , with X=Zn or S (see text) in the UVES GR B-DLA absorber sa mple (red diamonds). Similar measurements 
in the sample of QSO-DLAs observed with UVES (jNoterdaeme et elLI I2008T ) are shown with squares. Black squares are 
for H2-bearing QSO-DLAs and blue ones for H2 non-detections. There is no H2 detection in the GRB-DLA sample. 
Right panel: Same as before but as a function of depletion factor, [X/Fe]. Due to their very low H° column densities, 
GRB 060607 and GRB 080310 are not featured in these plots. 



a handle on the dust content of the systems, iron was se- 
lected as a proxy for refractory elements. Total Fe + column 
densities were calculated by summing up the contributions 
of both the ground state and the metastable energy levels 
of Fe + when detected. The contribution of the latter levels 
never exceeds 0.2 dex. In addition, we note that the cor- 
responding absorption lines are not detected at all in the 
DLAs towards GRB 050820 and GRB 080413A. 

H2 is not detected in any of the systems of the sample. In 
Table [31 we give 3a upper limits on the molecular hydrogen 
column densities and corresponding H2 molecular fractions, 
/ = 2N(E 2 )/(2N(rl 2 ) + N(R )). We also give 3ct upper 
limits on N(C°). 

Apart from the DLA towards GRB 050820 where the 
metallicity is as high as [Zn/H] = —0.4, the GRB host 
galaxies in the UVES sample have low metallicities, —2.2 < 
[X/H] < —1.3. In addition, the depletion factors are usually 
quite small, as expected for metal-poor environments hav- 
ing a low dust content. The only exception to this is the 
GRB 050922C-DLA for which sulphur is used as a proxy 
for non-refractory elements. In this particular case, the S II 
line profile does not follow that of Fe 11 at all, exhibiting an 
inversely asymmetric shape which peaks at the velocity of 
the detected Fen* absorption. Significant overabundance 
of a-clements and/or ionization effects might be present 
in this system and bias the observed [S/Fe] ratio high. 
More generally, the [S/H] and [S/Fe] ratios given in Table [3] 
should be considered cautiously as they are strictly speak- 
ing upper limits to the true metallicities and dust depletion 
factors. This makes the result that both metallicities and 
dust depletion factors in the sample of GRB-DLAs observed 
with UVES (and currently all other high-resolution spec- 
trographs) are usually modest, even stronger. Interestingly, 
this is most clearly seen in the three systems exhibiting the 



highest neutral hydrogen column densities, log JV(H°) ~ 22, 
where [X/H] < -1.5 and [X/Fe] < +0.2. The nature of 
these systems is discussed in Sect. 13.31 

The observed H° column density distribution of Swift- 
era, GRB-DLAs is relatively flat a nd not highly skewed 
towards extremely high values (see iJakobsson et al.l l2006. 
and also Table [3]). In addition, because there is mount- 
ing evidence that the bulk of the H° gas in GRB-DLAs 
has a galactic origin, with distances from the burst larger 
than 100 pc up to kiloparsecs (see Sects. 12.31 and 13. ip . 
the comparison of GRB-DLA properties to results drawn 
from QSO-DLA samples must be relevant. This is corrob- 
orated by our finding that the DLA cloud in the case of 
GRB 050730 is smooth and diffuse, and not apparently per- 
turbed by star-forming regions, with a broadening parame- 
tre and physical size that are typical of the Galactic ISM. 
In Fig. [8l we compare the H° column densities with metal- 
licities and dust depletion factors for the UVES GRB-DLA 
sample with QSO-DLAs. Studies of H 2 in large QSO-DLA 
samples observed with VLT /UVES ([Ledoux et all 120031 : 
iNoterdaeme et al"1l2008l ) have shown that the presence of 
H 2 does not strongly depend on the total H° column den- 
sity, and that the probability of finding H2 with high molec- 
ular fractions, log/ > —4.5, is only marginally higher for 
log7V(H°) > 20.8 than below this limit (19% and 7% re- 
spectively). The probability of H2 detection is of the or- 
der of 10% for the entire QSO-DLA population, but is 
actually strongly dependent on the metall icity, rising up 
to 35 % for systems with [X/ H] > —1.3 (jPetitjean et al.l 
2006; INoterdaeme et al.ll2008l ). In contrast, only about 4% 
of the [X/H] < —1.3 systems have log / > —4.5. Applying 
this number statistics to the UVES GRB-DLA sample 
(see Table [3J), composed of one [X/H] > —1.3 and four 
[X/H] < — 1.3 DLAs, leads to a binomial probability of de- 



tecting H2 at least once of 45%. This is consistent with the 
fact that no detection was achieved in practice. 

To estimate the dust content of DLAs, it is useful to 
calculate the column density of iron in dust, iV(Fe)dust = 
(1 - iQ-[x/Fc] )MXVFe/xl fly ladilo et al.ll2006h. In their 
survey of 77 QSO-DLAs, INoterdaeme et all <|2008h have 
shown that all H 2 -bearing systems have log iV(Fe)dust > 
14.7. In addition, about 17% of the systems with 14.7 < 
log A r (Fe)dust < 15.2, and 55% of the systems above this 
limit, have detected H2 lines with log/ > —4.5. The 
log 7V(Fe)dust values pertaining to GRB-DLAs are given in 
Table [3l Except in the case of the high-metallicity DLA 
towards GRB 050820, which exhibits logA(Fe) dus t ~ 16, 
GRB-DLA values lie either within (in three cases) or below 
(in one case) the above range. Therefore, from QSO-DLA 
statistics the binomial probability of not detecting H2 at 
all in the UVES GRB-DLA sample is 26%. From a 1.1a re- 
sult, relying on the existence of a single system, it would be 
premature to conclude that GRB-DLA observations are in- 
consistent with known Q SO-DLA propert i es. Th is is in con- 
trast with the claim bv iTumlinson et "all (|2007f L based on 
an even smaller and inhomogcncous sample of three high- 
and two medium-resolution GRB afterglow spectra, that 
H2 is generally deficient in GRB-DLAs compared to QSO 
lines of sight. 



3.3. Gas physical conditions 

The current lack of H2 detection along the lines of sight 
to GRBs might be related to the physical nature of DLAs 
exhibiting the highest neutral hydrogen column densities, 
logA(H°) ~ 22. Indeed, these systems make up about 
half of the UVES GRB-DLA sample (i.e., the systems 
towards GRB 050730, GRB 071031 and GRB080413A). 
Interesting ly, three systems in the la rge QSO-DLA sample 
studied by INoterdaeme et all (|2008h share similar proper- 
ties with no detected H2 but some of the highest column 
densities of iron in dust (log A(Fe)dust > 15.5) as well as the 
highest neutral hydrogen column densities ever measured 
along QSO lines of sight (i.e., log N(B°) = 21.70 towards 
Q 0458-0203, logA(H°) = 21.80 towards Q 1157+0128, 
and logA^H ) = 21.40 towards Q 1209+0919). The much 
smaller incidence rate of such DLAs in QSO samples 
is probably due to the narrow range of galactocentric 
radii probed by GRB afterglow observations. The prop- 
erties of these systems are at variance with Galactic 
lines of sight, where clouds with log A^(H°) > 21 usu- 
ally have 1ok7V(H 2 ) > 19 (see, e.g.. ISavage et all fl977t 
I Jenkins fe Shavalll979fl . 



ITumlinson et alJ (|2007l ) argue that UV radiation fields 
from recently formed hot stars with intensities 10-100 times 
the Galactic mean value must be invoked to explain GRB- 
DLA observations. This is conceivable if the clouds are em- 
bedded in or are l ocated close to re gions of intense star 
for mation (see also IChen et al.l [2007t) . However, as noted 
by ITumlinson et al.l ( 2007h this contrasts with the envi- 
ronment of massively star-forming regions like 30 Doradus 
in the LMC where the existence of clouds with high Hg 
molecular fractions i s ubiquitous (jBluhm fc de Boerll200ll : 
Dan forth et al.ll2002h . In other words, the presence of in- 
tense UV radiation fields does not necessarily prevent H 2 
from forming. 



In addition to a low dust content (see Sect. 13. 2| ). and the 
correspondingly low formation rate of H2 onto dust grains, 
another reason why H2 is not detected in those QSO- and 
GRB-DLAs with high H° column densities could be the 
particle density in the neutral gas being too small. For a 
cloud at equilibrium, by equating the H2 formation (i?dust) 
and photo-dissociation (i? diss) rates, one can write (see, e.g., 
Iffirashita fe Ferrarall2005T) : 



-Rdi ss ^V(H 2 ) = nn-RdustA^H) 



(4) 



where A(H) is the total hydrogen column density and 
A(H2) the molecular hydrogen column density, which then 
linearly depends on the particle density, hh ■ Therefore, the 
lower the density the lower the molecular fraction. In the 
case of the DLA towards GRB 050730, from its high H° 
column density, logA(H°) = 22.10, and considering the 
best-fit linear cloud size of I = 520±^° pc, a meaningful 
order of magnitude estimate of the particle density can be 
inferred: nu » 5 — 15 cm" 3 . Using E q. H] and the prescrip- 
tions bv iHirashita fc Ferraral ((2005!), we can estimate the 
UV intensity at hv = 12.87 eV outside the GRB 050730- 
DLA cloud, averaged over the solid angle, taking into ac- 
count H2 self-shielding. In order to keep the amount of H2 
below the observational threshold, logA(H2) < 13.8 (see 
Table \S§ , the UV intensity must be at least of the order of 
the Galactic radiation field: Jlw ^3x 10~ 20 erg s _1 cm~ 2 
Hz -1 sr _1 . This differs by more than a factor of ten from the 
results by Tumlinson et al. (see above). This is because in 
the latter study a scaling of the Galactic H2 formation rate 
based on the metallicity rather than on the dust-to-gas ra- 
tio, k = ioi x / H ](i-io[ Fc / x ]), was assumed or, equivalently, 
all metals are incorporated into dust grains. However, a 
strong correlation between depletion f actor and metallicit y 
is observed in QSO-DLAs (see, e.g.. iLedoux et all 12003). 
with smaller depletion factors in low metallicity systems, 
indicating that the dust-to-metal r atio increases markedly 
in the course of chemical evolution (|Vladildl2004l ). 

While, as shown above, particle densities are probably 
at most moderate in GRB-DLAs, with nn ~ 5 — 15 cm~ 3 , 
the ambient UV radiation field intensities needed to explain 
the lack of H2 in systems with high H° column densities are 
only of the order of 1-10 times those observed in the Galaxy, 
thus not particularly high. As shown in Sect. I3T21 the main 
factor driving the absence of H2 in the current sample of 
GRB-DLAs observed with VLT/UVES is a combination of 
low mctallicitics and low dust contents. 



4. Summary and prospects 

In this paper, we have presented a detailed analysis of 
the DLA cloud (logiV(H ) = 22.10) at the GRB 050730 
host-galaxy redshift of z a b s = 3.969 using high-quality 
VLT/UVES spectra. Accurate metal abundances are de- 
rived showing that this GRB-DLA has both low metallic- 
ity, [S/H] = —2.18 + 0.11, and no dust. In addition, we 
measured [N/S] = -0.98 + 0.05 and [Ar/S] = -0.11 + 0.06. 
The former ratio is indicative of a high nitrogen enrich- 
ment dominated by intermediate-mass stars (4 — 8 M@), 
before significant oxygen production b y supernovae, and a 
past l ow star-formation effi ciency (e.g. jHenry fc Prochaskal 
120071 : iPetitiean et al.l l2008). The latter ratio is a tra cer of 
the radiation field in H° regions (|Vladilo et al.ll2003D . The 
observed high abundance of argon, typical of the local ISM, 



is consistent with the expectation for a neutral medium em- 
bedded in a soft ionizing continuum. 

From the self-consistent photo-excitation modelling of 
absorption lines from electronic transitions to an unprece- 
dentedly large number of Fe + energy levels, the distance 
of the absorbing neutral material to GRB 050730 is con- 
strained to be d w 0.5 kpc. Similar photo-excitation mod- 
elling applied to the observed column densities and upper 
limits of S 3+ energy levels has led to a l ower limit on th e 
distance of the S 3+ -bearing gas of 0.4 kpc (|Fox et al.ll200l . 
From detailed calculations, we find that at such distances 
GRB afterglow photons cannot ionize H° nor dissociate 
H2 significantly. From the above Fe + column density mod- 
elling, a linear cloud size of I = 520l^gQ pc is inferred, sug- 
gesting that GRB-DLAs typically exhibiting very high H° 
column densities are diffuse metal-poor atomic clouds with 
large physical extents, I ^5 100 pc, and, at most, moderate 
particle densities, nn ~ 5 — 15 cm . 

In order to understand the lack of H2 further, we have 
built up the most up-to-date sample of Swift-eia. z > 1.8 
GRB-DLAs observed with VLT/UVES. We showed that 
the current non-detections of H2 in GRB-DLAs are con- 
sistent with QSO-DLA statistics where H2 is present with 
high molecular fractions, log / > —4.5, in only 10% of the 
global QSO-DLA population. Moreover, the lack of H2 in 
the GRB-DLA sample can be entirely explained by the low 
metallicities, [X/H] < —1, and low dust contents of the sys- 
tems. With the possible exception of the high-metallicity 
DLA towards GRB 050820, there is no need for enhanced 
UV radiation fields from recen tly formed hot stars to ex- 
plain ^_thi3_ljj£k_of IL2 (but see iTumlinson et all 120071 : see 
also IChen et al.l 120091 who derived high mean UV radi- 
ation field intensities for four bright GRB host galaxies, 
whose associated DLAs all have fairly high metallicities). 
The additional non-detection of Ci absorption lines, with 
log A(C°)/A(S + ) < —3, suggests a warm neutral medium, 
with kinetic temperatures T^ n <; 1000 K, as found in most 
QSO-DLAs (ISrianand et al.ll2005f ). 

The evidence for low to very low GRB-DLA metallici- 
ties coming from high spectral resolution observations is in 
contrast with the current sample of metallicity and deple- 
tion measurements drawn fro m low- and/or inte r mediate- 
resolution spectroscopy (e.g. , iFvnbo et al.1 120061 : ISavagliol 
l2006t iProchaska et al.ll2007bf ) which points towards higher 
metallicities than in the QSO-DLA population (see also 
IFvnbo et al.l I2008L who tried to reconcile the differ- 
ent metallicity distributions of QSO- and GRB-DLAs). 
Moreover, hints at the presence of H2 and/or Ci lines have 
been reported at th e host-galaxy r edshifts of GRB 50401 
([Watson et alll2006f). GR B 060 206 IFynbo et al.l l200l. and 



GRB 070802 (jEhasdottir et al.l l2009h. while the metallici- 
ties of these hosts are amongst the highest measured in 
GRB-DLAs. This suggests that it is only a matter of time 
before H2 and/or other molecules are detected beyond any 
doubt in GRB-DLAf(£|. Care should be exercized, however, 
not to overinterpret the current inhomogeneous body of 
data, as claims for high Zn + column densities and deple- 
tion factors of zinc compared to iron are based on a handful 



of absorbers at intermediate redsh ifts, 1 < z < 2, for which 
A(H°) measurements are lacking (|Savaglioll2006h . 

The VLT/UVES sample of high-redshift GRB after- 
glows shows a lack of GRB-DLAs with both high H° col- 
umn density and high metallicity (see Fig. [SI left panel). 
However, this does not mean that such systems do not 
exist and an anti-correlation between absorber dust con- 
tent and afterglow optical brightness should be checked. 
Indeed, the associated extinctio n will b e high in these 
systems (|Vladilo fc Perouxl 120051 : see also IProchaska ^it al.l 
20091 ) and their observation with high-resolution spectro- 
graphs, even with the advantage of fast response, possi- 
bly difficult. This is consistent with the findings of multi- 
colour photometry that a large fraction of burst s are dark 
because of high dust extinction in their hosts dRol et al 



20071: Uaunsen et all 120081 : iTanvir et~aTI 120081 : ICenko et al 



6 Since the submission of this paper, CO molecules with ex- 
ceptionally high column densities have been detected at the 
host-galaxy redshi ft of GRB 080607 in a low-resolution Keck af- 
terglow spectrum (IProchaska et al]|2009D . The absorbing gas is 
estimated to have roughly solar metallicity. 



2009). If a bias exists in the UVES sample, it could also 
be against GRB-DLAs located at much smaller distances 
from their bursts than currently found because, as we have 
shown, high metal (resp. H2) column densities are able to 
survive the ionizing (resp. dissociating) effects of the inci- 
dent GRB afterglow radiation. 

This potential problem should be alleviated by the 
second-generation VLT instrument X-Shooter, currently in 
its commissioning phase, to complement the capabilities of 
UVES at the VLT, Unit 2 - Kueyen, telescope. X-Shooter 
will be the instrument of choice to cover any redshift from 
z = 0.3 up to z — 10 for 1600 A (rest frame) Fen and other 
metal lines in GRB afterglows fainter, or redder, than ob- 
servable with UVES. In the meantime, we plan to apply the 
techniques presented in this paper for GRB 050730 to other 
GRB lines of sight observed with UVES to start building up 
a sample of burst/DLA distance determinations and map 
the distribution of neutral gas in GRB host galaxies. 
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